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ABSTRACT 

We present BVRI photometric and low-resolution spectroscopic investigation of a 
, type II core-collapse supernova (SN) 2008gz, which occurred in a star forming arm and 

within a half-light radius (solar metallicity region) of a nearby spiral galaxy NGC 3672. 
^sO \ The SN event was detected late, and a detailed investigation of its light curves and 

spectra spanning 200 days suggest that it is an event of type IIP similar to archetypal 
SNe 2004et and 1999em. However, in contrast to other events of its class, the SN 2008gz 
£S| ', exhibits rarely observed V magnitude drop of 1.5 over the period of a month during 

■ plateau to nebular phase. Using 0.21 mag of Ay as a lower limit and a distance of 25.5 

Mpc, we estimate synthesized 56 Ni mass of 0.05 ± 0.01M Q and a mid-plateau My of 
— 16.6 ± 0.2 mag. The photospheric velocity is observed to be higher than that was 
observed for SN 2004et at similar epochs, indicating explosion energy was comparable 
to or higher than SN 2004et. Similar trend was also seen for the expansion velocity of 
H-envelopes. By comparing its properties with other well studied events as well a s by 
J-j ' using a recent simulation of pre-SN models of iDessart. Livne fc Waldmanl (|2010T ). we 

infer an explosion energy range of 2 — 3 x 10 51 erg and this coupled with the observed 
width of the forbidden [O i] 6300-6364A line at 275 days after the explosion gives an 
upper limit for the main-sequence (non-rotating, solar metallicity) progenitor mass of 
17 Mq. Our narrow-band Ha observation, taken nearly 560 days after the explosion 
and the presence of an emission kink at zero velocity in the Doppler corrected spectra 
of SN indicate that the event took place in a low luminosity star forming H II region. 

Key words: supernovae: general — supernovae: individual: SN 2008gz — galaxies: 
individual: NGC 3672 
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1 INTRODUCTION 

Core-collapse supernovae occur in late type galaxies and 
their locations are usually associated with the regions of 
high stellar surface brightness or recent /ongoing star for- 
mation, suggesting that they re present the end stages of 
massive stars (M > 8-10 Mq) (|Anderson fc Jame 3 120091 : 

* e-mail: roy@aries.res.in, rupakroyl980@gmail.com 



lHakobvan et all |2009T ). Observationally, these events are 
classified into H-rich type II SNe which show prominent 
H- lines in their optical spectra, and H-deficient type Ib/c 
SNe which don't show the trace of H-lines. Ic events lack 
He-lines as well. Type II SNe constitute about 70% of all 
the c ore-collapse SNe (|Cappellaro et al J 1 19991 : | Smith et al.l 
2009) and their light curves and spectra differ significantly. 
In type IIP, the optical light remains constant for about 
100 days (called the plateau phase) and then decays expo- 
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nentially. The spectra are featured due to strong P-cygni 
profiles, while in type IIL SNe, a linear decline in its optical 
light and strong emission lin es are observed. Type Hn events 
show narrow emission lines |Filippenkol 1 19971 ; ISmartt et al.l 
l2009al ). 

Theoretically, the explosion mechanism consists of the 
collapse of the progenitor star's Fe-core, formation of a shock 
wave, the ejection of stellar envelope and formation of a neu- 
tron star or a black hole. The shock wave generated through 
the reversal of core-collapse, breaks out the stellar surface 
of progenitor as a hot fireball flashes in X-ray and ultravi- 
olet continuing from a few seconds to a few days. In H-rich 
events, the shock-heated expanding stellar envelope cools 
down, by recombination of H and sustains the plateau phase 
of IIP SNe, while the post-maxima/plateau light curves are 
powered by the radioactive decay of 56 Co into 56 Fe. Though 
the explosion mechanism is similar to these events, they dif- 
fer largely in energetics, e.g. IIP SNe are observed to form 
a sequence from low-luminosity, low- vel ocity, Ni-poor events 
to bright, high- velocity, Ni-rich objects (|Hamu v 2003). Thus, 
a detailed investigation of individual core-collapse events is 
important for understanding the nature and environment of 
progenitors. They generally probe the star formation pro- 
cesses, galactic chemical evolution and constrain the stel- 
lar evolutionary models jHeger et al. 120031 ; ISmarttl l2009bl ; 
lHabergham. Anderson k, James! 120101 ). Type IIP SNe also 
turned out to be good standardizable cand les (|Hamuvll2002l ; 
IPoznanski et al.ll2009l ; lOlivares etailbOld ) . 

The SN 2008gz event was discovered on November 5.83 
UT, 2008 by Koichi Itagaki using a 0.6m telescope in the spi- 
ral galaxy NGC 3672 at an unfiltered magnitude of about 
16.2. On November 7.84 and 8.84 UT, an independent dis- 
covery of this new transient was reported by R. Martin 
from Perth Observatory as a part of "Perth Automated 
Supernova Search Program" by using 0.6m Lowell Tele- 
scope. The red magnitude of this new object was about 
15.5 dNakano fc Martini 120081 ). On November 11.25 UT, 
iBenetti et alj l|2008h took the first spectra of this event with 
the 3.5m TNG (+ DOLORES; range 340-800 nm, resolu- 
tion 1.0 nm) and showed tha t it is a type II supernova event 
and by using GELATO tool (|Harutvunvan et al.ll2008h , they 
found that the spectrum of SN 2008gz resembles best with 
that of a H-peculiar event SN 199 8A, taken at 62 days after 
explosion (Pas torello et al.ll2005l ). Assuming the recession 
velocity of the host galaxy ~ 1862 kms -1 , they found the 
expansion velocity of hydrogen layer was about 6600 kms -1 . 
An independent regular BVRI CCD photometric monitor- 
ing of SN 2008gz was carried out since November 10, 2008 
by using lm Sampurnanand Telescope at Nainital, India. 
We also collected spectra in optical (0.4 — 0.9 /im) with 2m 
IUCAA, India; 3.6m NTT, Chile; 6m BTA, Russia; 3.6m 
TNG, Spain. 

In this work, we present results of optical photometric 
and low-resolution spectroscopic investigation of SN 2008gz. 
We adopt time of explosion to be August 20.0, 2009 or JD 
2454694.0 having uncertainty of a few days (see ^2.2l for de- 
tails). Hence the time of post/pre-explosion are rounded off 
to nearest day and they are referred with + and — signs 
respectively. Basic properties of SN 2008gz and its host 
galaxy NGC 3672 are given in Table [T] The paper is or- 
ganized as follows. Sj2] and fj3] present the photometric and 
spectroscopic observations and a brief description of light 



Table 1. Properties of the host galaxy NGC 3672 and SN 2008gz. 



Parameters 



Value 



Ref. a 



NGC 3672: 

Type 

RA (J2000) 
DEC (J2000) 
Abs. Magnitude 

Distance 
Scale 

Distance modulus 

Apparent radius 
Inclination angle 
Position angle 
Heliocentric Velocity 

SN 2008gz: 

RA (J2000) 
DEC (J2000) 

Location 

Deprojected radius 
Explosion epoch (UT) 
Discovery date (UT) 



SAc 

a = ll h 25 m 2?48 
6 = -09°47'43'.'0 
Mg = -20.59 mag 

D = 25.5 ± 2.4 Mpc 

1" ~ 123 pc, 1' ~ 7.4 kpc 

H = 32.03 ±0.21 

r 25 = l.'4(~ 10.4 kpc) 
inc = 56.2° 



On 



: 6.5° 



c^helio = 1864 ± 19 kms -1 

a = ll h 25 m 3?24 
5 = -09°47'51'.'0 

13" E, 7" S 

r SN = 23'.'37 (~ 2.81 kpc) 

20.0 August 2008 
(JD 2454694.0) 
5.83 November 2008 



1 
1 
1 
1 

m 



a (1) Hyper LEDA - http://leda.uniy -lyonl.fr; (2) 
iNakano fc Martini feOOSl) 

curves and spectra. In ij4]we study evolution of some impor- 
tant line profiles, where as in SJH the velocity of the photo- 
sphere and the H-ejecta is estimated by using SYNOW code 
jBranch. Baron fe Jeffery||200ll ; iBranch et atll2002l ; iBaron I 
2005) that describes spectroscopic observations. Distance, 
extinction and the evolution of colour and bolometric lumi- 
nosity are studied in ijfj]and ^respectively. The amount of 
synthesized 56 Ni mass, environment and energetics of the 
progenitor are estimated and discussed in jj8] We also made 
a comparative study of this event with other type IIP SNe 
in ij9] At last a summary is presented in i]10l 



2 BROAD BAND PHOTOMETRY 
2.1 BVRI data 

Initial pre-SN images (-403d) of the host galaxy NGC 3672 
at VRI bands were obtained from Perth Observatory, as a 
part of supernova search program for another type la SN 
2007bm, which occured in the same galaxy. Images were 
taken with 512 x 512 CCD camera mounted on a 0.6m Low- 
ell Telescope, covering around 5x5 square arcmin on the 
sky. The FWHM seeing was about 2'.' 5. 

SN 2008gz was observed at different epochs from differ- 
ent observatories around the world. The major part of moni- 
toring was carried out in Johnson BV and Cousins RI bands 
from lm Sampurnanand Telescope (STj3 at the Aryab- 



We used a 2048 X 2048 CCD camera having a square pixel 
of 24 ^m mounted at the f/13 Cassegrain focus of the telescope. 
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Figure 1. SN 2008gz in NGC 3672. V band image from 1 m ST, 
India. Area of about 8x8 arcmin is shown, with location of SN 
marked with a cross and reference standard stars marked with 
circles. 



hatta Research Institute of Observational Sciences (ARIES), 
Nainital, India. SN 2008gz was observed during 10 Novem- 
ber 2008 (+87d) to 17 May 2009 (+275d). We could not 
detect SN in the observations of 19 November (+462d) in 
B, of 13 February 2010 (+547d) in VRI and of 14 February 
2010 (+548d) in BVRI. In addition to lm ST, observations 
of SN 2008gz at BVRI bands were also obtained on 21, 22, 
24 and 25 March 2009, with IFOSC mounted on 2m IGO, 
IUCAA, India and on 17 May 2009 with EFOSC2 mounted 
on 3.6m NTT, ESO, Chile. The journal of observations is 
given in Table [ffl 

Photometric observations included acquisition of sev- 
eral exposures with exposure time varying from 100 to 300 s 
in different filters. Several bias and twilight flat frames were 
obtained for the CCD images. Bias subtraction, flat fielding, 
cosmic ray removal, alignment and determination of mean 
FWHM and ellipticity in all the object frames were done us- 
ing the standard tasks available i n the data reduction soft- 
wares IRAF^ and DAOPHOf^ ( Stetson! [l987l , Il992l ). The 
FWHM seeing at V band varied from 2" to 4", with a me- 



Plate scale of the CCD chip is 0.38 arcsec per pixel, and the entire 
chip covers a field of 13 X 13 square arcmin on the sky. The 
gain and readout noise of the CCD camera are 10 e~ per analog- 
to-digital unit and 5.3 e~ respectively. All the observations were 
carried out in the binning mode of 2x2 pixel. 

2 Table [2] is only available in electronic form. 

3 IRAF stands for Image Reduction and Analysis Facility dis- 
tributed by the National Optical Astronomy Observatories which 
is operated by the Association of Universities for research in As- 
tronomy, Inc. under co-operative agreement with the National 
Science Foundation. 

4 DAOPHOT stands for Dominion Astrophysical Observatory 
Photometry. 
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dian value of around 2'.'5. About 10% of the images taken 
at large zenith distance had highly elongated PSF (ellip- 
ticity > 0.2). For final photometry, we co-added the indi- 
vidual frames to increase the signal-to-noise ratio. The pre- 
processing steps for images taken from other than lm ST 
were also performed in similar fashion. 

Fig. [T] shows the location of SN 2008gz in the galaxy 
NGC 3672. The SN flux is expected to have substantial con- 
tribution from the host galaxy background, due to its prox- 
imity to the galaxy centre, its location in a spiral arm and 
a high inclination angle (56.2°; Table [1} of the galaxy. At 
early phases, SN flux dominates the total flux, thus with a 
PSF-fitting method we were able to remove the galaxy con- 
tribution. At later epochs (e.g. end of plateau in case of type 
IIP events), galaxy flux may brighten the SN light curves 
by 0.5 to 1 mag depe nding on its locati on in the galaxy 
llPastorello et all 120051 ). We used ISlffi (|Alard fc Luptonl 
Il998h to get galaxy template subtracted flux of supernova. 
As a template we used the BVRI images taken on 14 Febru- 
ary 2010 (+548d) from lm ST, India in good seeing con- 
ditions. We note that the galaxy subtraction using pre-SN 
(-403d) VRI images taken from Perth Observatory, gave no 
SN contribution above noise level in our images recorded on 
19 November 2009, 13 and 14 February 2010. In order to 
verify ISIS results, we also performed the galaxy template 
subtraction scheme independently using self-written scripts 
employing IRAF tasks which included alignment, PSF and 
intensity matching of the galaxy template and SN images, 
and subtraction of template from SN images. Fig. [2] shows 
images with and without template subtraction. PSF-fitting 
method was applied on the subtracted images. Our magni- 
tude was found to be consistent with the ISIS ones having a 
typical scatter of ~0.1 mag in BVRI; this scatter is of the 
order of 0.05 mag in plateau phase (see Fig. [3?h. 

In order to cali brate ins t rumen tal magnitudes of SN 
2008gz, we observed lLandoltl (|2009l ') standard fields SA92 
and PG0231 in BVRI with lm ST on 15 November 2008 
under good night conditions (transparent sky, FWHM see- 
ing in V ~ 2"). The data reduction of SN and Landolt 
fields were done using profile fitting technique and the in- 
strumental magnitudes were converted into standard system 
fo llowing l e ast-sq uare linear regression procedures outlined 
in IStetsonl <lT992h . We used mean values of atmospheric ex- 
tinction coefficients of the site viz. 0.28, 0.17, 0.11 and 0.07 
mag p er unit airmass for the B, V , R and / bands respec- 
tively l|Kumar et al.1 2000). A set of 13 stars having a colour 
range of —0.33 < (B — V) < 1.45 and brightness range of 
12.77 < V < 16.11 were used to derive the following zero 
points and colour coefficients: 

b — B + (5.37 ± 0.02) + (0.01 ± 0.02)(B - V) 

v = V + (4.89 ± 0.01) + (-0.02 ± 0.02)(B - V") 

r = R+ (4.72 ± 0.01) + (0.02 ± 0.02)(V - R) 

i = I+ (5.17 ± 0.02) + (0.02 ± 0.02)(V - I) 

Here B, V, R, I are the standard magnitudes and b, v, r, 
i are corresponding instrumental magnitudes corrected for 

5 http: / /www2. iap.fr/users/alard/package. html 

6 Figure [3] is only available in electronic form. 
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Table 2. Journal of photometric observation of SN 2008gz. The full table is available online. Please see the supporting information 
section for detail. 



UT Date 


JD 


Phase 


B 


V 


R 


I 


Telescope 


Seeing 


Ellipticity 


(yy/mm/dd) 


2454000+ 


(day) 


(s) 


w 


(s) 


(s) 




(") 




2007/07/09 


290.94 


-403 




240 


180 


180 


LT 


2.5 


0.05 


2008/11/10 


781.48 


+87 


2x300 


2x300 


2x250 


2x250 


ST 


2.4 


0.18 
























Figure 2. Template subtraction for SN 2008gz. V-band image from 10 December 2008 (+ 117 d) with SN is shown in leftmost panel, 
the template frame from 14 February 2010 (+ 548 d) without SN is shown in middle, while the rightmost panel shows the template 
subtracted image. All the images are around 3 X 4.5 arcmin. North is up and East is to left. 



time and aperture. A typical scatter in the photometric so- 
lutions to the Landolt standard stars is found to be ~ 0.03 
mag for BVRI. Table lists the calibrated magnitudes for 
a set of ten stable secondary standards in the SN field, while 
calibrated _BV".R7magnitudes of SN 2008gz are presented in 
Table 3] For SN, we quote ISIS derived errors (lcr uncer- 
tainty), which is consistent with the RMS scatter in the 
magnitude of standard stars determined from night-to-night 
repeatability over entire period (~ 215d) of SN monitoring. 
Large errors in 2m IGO and 3.6m NTT data arises due to 
mismatch in the PSF and pixel scale. 

2.2 Optical light curve 

Fig. H] shows BVRI light curves of SN 2008gz ranging 
from +87d to +275d since the time of explosion. We 
also present the light curves of other we ll studied nearby 
(D < 12Mpc) typ e IIP SNe yiz. 2004et ilSahu et al 1 l2006k 
iMisra et all 120071) . 1999em l|Elmhamdi et all 120033 ) and 
1999gi (|Leonard et al.ir2002th scaled in time and magnitude 
to match the transition between plateau to nebular phases. 
It is seen that SN 2008gz was detected close to the end of 
its plateau phase and its light curve resembles well with 
the above three template IIP events and hence we could 
determine (and adopt) the time of inflection (plateau to 
nebular), ti of 115 ± 5 d, by adjusting the template light 
curves to get the best match to SN 2008gz data points. 
This derived plateau dur ation is typical for type II P events 
|Elmhamdi et al.ll2003bh and it is also consistent with the 
fact that SN 2008gz was not visible around three months 
before the discovery date (Novem ber 5.83, 2008) at the leve l 
of unfiltered magnitude of 19.0 (jNakano fc Martini [2008). 
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SN 1999em 
- SN 1999gi 
SN 2004et 
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Figure 4. Light curve in BVRI magnitudes of SN 2008gz. The 
light curves are shifted for clarity, while for other SNe, they are 
scaled in magnitude and time to match with SN 2008gz. 



Further, similarities between bolometric tail luminosity of 
SN 2008gz (see with that of SN 2004et and SN 1999em 
indicates that probably explosion happened about 82 days 
before the discovery date. Analysis of first spectrum of SN 
2008gz (see £|3.2[1 also revealed few similarities between the 
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Table 3. Identification number (ID), coordinates (a, 5) and calibrated magnitudes of stable secondary standard stars in the field of SN 
2008gz. Errors in magnitude represent RMS scatter in the night-to-night repeatability over entire period of SN monitoring. 



Star QJ2000 <5j2ooo B V R I 

ID (h m s) (° ' ") (mag) (mag) (mag) (mag) 



1 


11 


25 


11.27 


-09 


48 04.4 


15.36±0.01 


14.36±0.03 


13.99±0.01 


13.49±0.01 


2 


11 


25 


18.34 


-09 


47 03.0 


15.70±0.01 


14.92±0.02 


14.68±0.01 


14.25±0.01 


3 


11 


24 


55.62 


-09 


46 07.9 


16.50±0.01 


15.82±0.03 


15.64±0.01 


15.26±0.01 


4 


11 


25 


14.22 


-09 


44 04.7 


16.91±0.01 


16.15±0.01 


15.93±0.01 


15.50±0.01 


5 


11 


25 


10.47 


-09 


49 17.0 


17.00±0.02 


16.23±0.03 


15.98±0.01 


15.52±0.04 


6 


11 


24 


53.17 


-09 


47 25.1 


17.09±0.02 


16.42±0.03 


16.22±0.01 


15.83±0.03 


7 


11 


24 


56.67 


-09 


41 35.1 


17.80±0.02 


16.99±0.02 


16.66±0.01 


16.19±0.04 


8 


11 


25 


14.78 


-09 


44 54.3 


18.50±0.07 


17.70±0.07 


17.49±0.06 


17.12±0.03 


9 


11 


25 


11.40 


-09 


44 54.5 


18.28±0.04 


17.75±0.03 


17.62±0.02 


17.22±0.04 


10 


11 


25 


09.31 


-09 


47 16.4 


18.83±0.09 


17.92±0.07 


17.63±0.03 


17.08±0.06 



kinematical properties of its ejecta to that of SN 2004et ob- 
served nearly 80 days after the burst. We therefore adopt 
time of SN explosion to be 82 ± 5 days before the discovery 
date and this corresponds to burst time, to of JD 2454694.0, 
however, we note that based on the first spectrum (Novem- 
ber 11.25 UT, 2008) and its simil arity with +62d spectrum 
of SN 1998 A ijBenetti et al.ll200St ) suggesting a time of ex- 
plosion of nearly 56 days before the discovery date (corre- 
sponding to plateau phase of ~ 90d) cannot be ruled out. 

In the late plateau phase (~ +90d) flatness behaviour 
in RI and decline trend in BV are clearly seen, which are 
similar to other IIP events. The V magnitude drop of 1.5 
mag from plateau phase (V ~ 16 mag at +100d) to neb- 
ular phase (17.5 mag at +130d), is slightly lower than 2 - 
3 mag drop for a typical IIP event (|Qlivares et aljfeoiol ). 
This shallow decline which is seen in BRI also, indicat- 
ing production of large 56 Ni mass (see £|8. 1[) . In contrary 
to this, very steep brightness decline at V has also been ob - 
served, e.g. 4.5 mag for SN 2007od (|Andrews et al.ll2010D . 
The nebular phase starts at ~ +140d, and it roughly fol- 
lows the decay slope of 56 Co to 56 Fe: 0.98 mag (lOOd)" 1 . A 
linear fit to the tail from +150d to +275d gives the follow- 
ing decline rates [in mag (100d) -1 ]: jb ~ 0.51, -yv ~ 0.98, 
7_r ~ 1.12, 7/ ~ 1.13 at B, V, R, I which is typical to 
the values found for IIP SNe. The flattening seen in B band 
light curve, though non-conclusive due to large scatter of the 
measure ments, has also been obs erved in othe r events, e.g. in 
1999e m (lElmhamdi et alj|2003ah and 1987A (|Suntzeff et all 
I1988T ) until +400d. 



3 LOW RESOLUTION SPECTROSCOPY 
3.1 Data 

Long-slit low resolution spectra (~ 6 to 14 A) in the optical 
range (0.33 - 1.0 (im) were collected at eight epochs during 
+87d to +275d; five epochs from 2m IGO, and one epoch 
each from 3.5m TNG, 6m BTA and 3.6m NTT. Journal of 
spectroscopic observations are given in Table [5] 

At 2m IGO, observations were carried out using 
IFOSC (IUCAA Faint Object Spectrograph and Camera) 
mounted at the cassegrain end off/10 refl ector l| Gupta et al.l 
120021 ; IChakrabortv. Das fc Tandonll2005l ). Slit spectra were 
recorded using 2048 x 2048 EEV CCD camera with 13.5 
fim pixel, having a gain of 1.8 e~per analog-to-digital unit, 



and readout noise of 6.3 e~ . Grism 7 with peak sensitivity 
at 500 nm and a slit width of l'/5 were used. Calibration 
frames (bias, flats, arcs) and spectrophotometric flux stan- 
dards were observed on each night. For SN, usually slits 
were placed across the spiral arm so as to make proper sky 
background and in one case at +170d, the galaxy centre was 
also observed. Spectroscopic data reduction was done under 
IRAF environment. Bias and flatfielding were performed on 
each frames. Cosmic ray rejection on each frame was done 
by using Laplacial kernel detection (Ivan Dokkum i 2001). Im- 
ages were coadded to improve the signal-to-noise ratio and 
one-dimensional spectra were extracted using apall task in 
IRAF which is based on optimal extraction algorithm by 
lHornell|l986h . Wavelength calibration were performed by us- 
ing identify task and about 15-18 emission lines of He and 
Ar which were used to find a dispersion solution. Fifth order 
fits were used to achieve a typical RMS uncertainty of O.lA. 
The position of OI emission skyline at 5577A was used to 
check the wavelength calibration and deviations were found 
between 0.5 to lA and this was corrected by linear shift in 
dispersion. The instrumental FWHM resolution of 2m IGO 
spectra as measured from O 1 5577A emission skyline was 
found to lie between 6A to 10A (~ 322 - 510 kms^ 1 ). Flux 
calibration was don e using spectrophotometric fluxes from 
lHamuv et alj (| 19941 ) and assuming a mean extinction for the 
site. Synthetic magnitudes were estimated using spectra to 
verify the accuracy of flux calibration and it was found to 
be accurate within 0.05 mag. 

Spectroscopic data reduction for DOLORES on 3.6m 
TNG, EFOSC2 on 3.6m NTT, and SCORPIO on 6m BTA 
were done in similar fashion and at around 6000A it had a 
resolution of 10A, 14A and 12A respectively. 

3.2 Optical spectra 

Fig. [5] and Fig. [6] show the rest frame spectra of SN 2008gz , 
corrected for recession velocity (1864 kms -1 ) of the host 
galaxy NGC 3672. We identify all the spectral features as 
per previously published line identifica tions for IIP events 
^Leonard et al.M2002l ; ISahu et alJkood ). In Fig. E] end of 
the plateau phase (+115d) and beginning of nebular phase 
(+140d) is clearly evident in the spectral evolution. 

The late plateau phase (+87d and +115d) spectra are 
marked by strong P-Cygni features of Ha, 1 7700A, Na I D 
5890, 5896A, and singly ionised Sc, Ba, Ti, Fe atoms, while 
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Table 4. Photometric evolution of SN 2008gz. Errors in magnitude are derived from ISIS and it denote lcr uncertainty. 



UT Date 


JD 


Phase" 


B 




V 


R 






Telescope 6 


Seeing c 


(yy/mm/dd) 


2454000+ 


(day) 


(mag) 




(mag) 


(mag) 




(mag) 




(") 


9008 /ii /1 o 


781.47 


4-87 
TO I 


i 7 onci-|-n 


016 


1 K 071 _(-n 01 9 


1 ^ 47^1+0 
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" with reference to the explosion epoch JD 2454694.0 
6 ST : 1 m Sampurnanand Telescope, ARIES, India; IGO : 2 m IUCAA Girawali Observatory, IUCAA, India; NTT : 3.6 m New 

Technology Telescope, ESO, Chile 
c FWHM of the stellar PSF at V band 
d Flat field problem 



the +140d and later spectra show significant drop in the ab- 
sorption strength of P-Cygni features. The spectra during 
+ 115d to +222d show the spectral evolution of the event 
from early to mid-nebular stage, while the last two spec- 
tra (+231d and +275d) are those typical shown during the 
late stages of a typical SNIIP. In Fig. [6l the +87d spectrum 
shows various atomic absorption lines over the weak con- 
tinuum. These lines are mainly due to elements present in 
the SN ejecta along with some earth atmospheric molecular 
lines (marked with ©) and absorption due NalD of Milky 
Way and host galaxy. On the other hand, the +275d spec- 
trum shows a typical nebular phase spectrum dominated by 
emission lines. 



Temporal evolution of P-Cygni nature of Ha is clearly 
seen viz. the emission component becomes narrower with a 
decrease in depth of associated absorption component dur- 
ing the transition of SN from plateau to nebular phase. 
The FWHM of emission component of Ha decreases from 
~ 5477kms~ 1 at +87d to ~ 3526kms~ 1 at +275d, indicat- 
ing decrease in opacity and temperature of Hi line emit- 
ting regions. For H/3, H7 and Hi5, the emission components 
are crowded with numerous metal lines. In Fig. [5] we also 
see impression of an additional P-Cygni component in the 
absorption profile of Ha and H/3. This is speculated as a 
footprint of high velocity emitting shells in SN ejecta. Sim- 
ilar signatures were also noticed in type IIP SNe 1999em 
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Table 5. Journal of spectroscopic observations of SN 2008gz. 



UT Date JD Phase" Range 

(yy/mm/dd/hh.hh) 2454000+ (days) /an 



Telescope 



Grating 
(gr mm -1 ) 



Slit width 
(") 



Dispersion 
(Apix- 1 ) 



Exposure S/N c 
(s) (pix- 1 ) 



2008/11/11/06.049 


781.88 


+87 


0.34-0.80 
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500 
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80 


2008/12/08/23.346 


809.46 


+ 115 


0.38-0.68 


IGO 


600 


1.5 


1.4 
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834.52 
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0.38-0.68 
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865.40 


+ 171 


0.38-0.68 
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1.5 
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913.27 
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0.38-0.68 


IGO 
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1.5 


1.4 


1800 
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916.33 
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0.38-0.68 


IGO 


600 


1.5 


1.4 
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7 


2009/04/03/22.643 


924.50 


+231 


0.61-1.00 
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550 


2.1 


3.5 


3x900 


70 
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969.04 


+275 


0.33-0.80 


NTT 


300 


1.0 


i.O 


2700 


16 


2009/05/17/01.308 


969.07 


+275 


0.55-1.05 


NTT 


300 


1.0 


4.2 


2700 


24 



a With reference to the burst time JD 2454694.0 
6 TNG : DOLORES on 3.5m Telescopio Nazionale Galileo, Italy; IGO : IFOSC on 2m IUCAA Girawali Observatory, India; BTA : 
SCORPIO on 6m Big Telescope Alt-azimuthal, Special Astrophysical Observatory, Russia; NTT : EFOSC2 on 3.6m New Technology 

Telescope, ESO, Chile. 
c At 0.6 ^tm 
d Only center of galaxy observed 



15 - 




+87 d 
+115d 
+140 d 
+171 d 




+218 d 
+222 d 

+231 d 
+275 d 



4000 6000 8000 

Rest Wavelength (A) 

Figure 5. Doppler corrected flux spectra of SN 2008gz from late 
plateau (+87 d) to nebular phase (+275 d). Prominent Hydrogen 
and metal lines are marked. 



and 2004et ijLeonard et al.ll2002l : lEi et al.ll2005l : ISahu et all 
I2006T ). 

The spectrum, labeled with +275d shows a typical 
late nebular phase spectrum marked by emission dominated 
permitted lines of Can 8498, 8542, 8662A and NalD as 
well as with the appearance of forbidden emission lines, 
i.e., [Oi]6300, 6364A; [Fen]7155A and [Can] 7291, 7324A . 
These forbidden lines are not observed at Earth due to high 
density of gas. [Call] is already seen in +87d spectrum, 
[Ol] appears at +171d, while the [Fell] appears in +275d 
spectrum. Increasing strength of these forbidden lines in- 
dicate expansion and rarefaction of SN ejecta with time. 



The [Fen] line is also noticed bv lPastorello et all (|2005t ) in 
+344d spectrum of SN 1998A and +346d spectrum of SN 
1987A, whereas for low lum inosity SN 1997D it was visible at 
+417d ()Benetti et al.l[200ll ). We determine relative strength 
/([Caii])/J([Feii]) of 2.06 for SN 2008gz, whereas for SNe 
1999em and 1987A, values of this ratio are respective ly 6.98 
and 24.68 at around +400d (jElmhamdi et al.fl2003ah . This 
indicates that the physical conditions of [Fen] formation in 
SN 2008gz may be similar to that of SN 1999em rather than 
SN 1987A. 



4 TEMPORAL EVOLUTION OF SPECTRAL 
LINES 

To further illustrate the nature of burst, we show in Fig. [7] 
(and describe below) the velocity profiles of Ha, H/3, NalD, 
Ban6142A and [O i] 6300, 6364A. The absorption dips and 
emission peaks due to SN are marked by downward and up- 
ward arrows respectively. For a spherically symmetric burst, 
the emission peak of P-Cygni profiles should be located at 
the rest wavelength of the corresponding line, while absorp- 
tion dip will be blueshifted reflecting the instantaneous ve- 
locity of corresponding line emitting region. In our rest frame 
spectra, Ha emission peak is found to be slightly blue shifted 
by ~ 406 kms -1 at +87d. Such blue shift in Ha emis- 
sion peak at the early e pochs was also observed f or other 
type II SNe (i.e., 19 87A jHamischik fc Dachslll987l) 1988A 
Turatto et all 19931). 19 90K dCappellaro et al.ll 19951 ) ,1993 J 



Matheson et al l l2000h. 



1999e m (lElmhamdi et al"1]2 003a) 



20051) 



1998A ilPastorello et all 

2005cs ijPastorello et all 
20091 1 . 2007od (|Andrews et al.ll2010t )). It is observed for few 
low lum inosity type II as well, e.g., SN 1999br in their early 
phases (IPastore lo et al.l|2004j) . On the basis of SN 1987A 



Chugail 1 1988) explained this phenomenon 



velocity profile, 
as an effect of diffused reflection of resonance radiation by 
the expanding photosphere. For SN 1987A at the +85d this 
velocity was about a few hundred kms _1 and hence compa- 
rable with SN 2008gz. On the other hand for SN 1998A and 
SN 1993J Ha emission peak velocity at comparable epoch 
were at least one order of magnitude higher than SN 2008gz. 
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Figure 6. Line identification for late plateau and deep nebular phases. 



At later epoch (+275d) Ha emission peak velocity for SN 
2008gz is about — 200 km s^ 1 which is considerably different 
from SN 1987A, where redshifted velocity of +600 km s" 1 
was observed. 

In most of our Ha profile, a narrow emission line 
(marked by 'D' in Fig. 0) is seen at zero velocity, which 
is probably due to underlying H n region and this is consis- 
tent with presence of low-luminosity H n region at SN loca- 
tion revealed by our Ha narrow band observation f H8.20 . 
Such a fea ture is not seen in well ob served nearby IIP 
SNe 2004et (|Sahu et alj2006h and 2004A (jGurugubelli et alj 
2008), which occurred in the outskirts of their host galaxies. 

Na I D 5892A P-Cygni nature is prominent in the spec- 
tra at all epochs with emission peak located at zero veloc- 
ity, indicating almost spherical distribution of Na ion in the 
ejected material. In high S/N spectra (+87d and +115d), 
two peculiar absorption dip in the SN NalD emission profile 
is seen at — 80kms -1 (marked by B), and at — 1630 kms" 1 
(marked by C), respectively. These are identified as NalD 
absorption respectively due to the interstellar matter in the 



host galaxy and the Milky Way (31} . In this regard, we note 
that due to highly inclined host galaxy (0i nc = 56.2°), true 
recession velocity for SN will be different from the adopted 
ones (1864 kms -1 ), and as measured from NalD absorption 
at SN location, it should differ by ~ 100 km s -1 , however, 
considering, low dispersion spectra, this will not change any 
of our conclusions. 

The absorption dips due to an s- process element Ball 
6142 A and for the element Sell 6248A are clearly visible in 
+87d and +115d spectra, however, they disappeared in the 
+ 140d spectrum. This is similar to luminous type IIP SNe 
2004et (|Sahu et al.| [2006) where these lines were prominent 
at +113d and barely ob servable at +163d and for 1999em 
l|Elmhamdi et al]|2003al ) it was visible up to +166d. In con- 
trar y to this, low luminosity IIP SNe 1997D , 2005cs and oth- 
ers (ITuratto et al.ll 19981 ; iBenetti et ai]|200ll ; iPastorello et all 
12004 ). these features sustain comparatively for longer time 
and observable till +208d. Low luminous type IIP SNe ex- 
pand with velocity, much slower than that of normal type 
IIP events. So, the Ba lines in low luminosity events sus- 
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4861 



-115d 




-275d 



Velocity (1000 km s" 1 ) 

Figure 7. Evolution of some important spectral lines of SN 2008gz, during its transition from late plateau to nebular phase. The dotted 
line at zero velocity corresponds to the rest wavelength. Upward arrows show emission peaks, while downward arrows show absorption 
dips. The feature 'A' is a high velocity P-Cygni component of Ha , features 'B' and 'C are due to NalD absorption from the Milky Way 
and the host respectively, and the emission feature 'D' is Ha emission due to underlying Hn region at SN location. The numbering 
sequences and resolution arc : 1 (+87 d/loA), 2 (+115 d/6A) 3 (+140d/12A), 4 (+171d/6A), 5 (+218d/10A), 6 (+222d/10A), 7 
(+231d/10A) and 7 (+275 d/14A). 



tain for longer time just because the ejecta takes more time 
to cool-down. For luminous type II-peculiar SNe 1987A and 
1998A, these feature are seen even at later epochs beyond 
+300d. For SN 2008gz, probably Ba abundance was low and 
s-process was not so effective like low luminosity events. 

Metastable [Ol] 6300, 6364 A lines start to appear at 
+ 171d, and its strength increases progressively in spectra 
at later epochs. A normalized profile for +275d is shown in 
Fig. \$\ which is nearly symmetric, indicating a spherically 
symmetric Oxygen ejecta. The average FWHM for [O i] lines 
is ~ 2500 kms -1 . Two component Gaussian fit results in 



Figure [8] is only available in electronic form. 



ratio of 7(6300)//(6364) ~ 1.8, which is at deviation from 
their strength ratio of 3 expected from transitional prob- 
ability for a rarefied gas at certain temperature. Smaller 
ratio for SN 2008gz may indicate higher opacity for 6300 
A in comparison to 6364A. We however note that the ra- 
tio of J(6300)//(63 64) is not always 3 for type IIP events. 
For the SN 1988A, ISpvromilid fl991) showed that at ini- 
tial epochs when optical depth of the ejecta is very high, 
value of / (6300)// (6364) ~ 0.952, whereas at late phase, for 
optically thin ejecta this ratio approaches to 3.03. 

Blueshift in [O i] lines, particularly at epochs later than 
+400d (i.e. in SNe 1999em and 1987A), is interpreted as an 
effective indicator of dust formation in the SN ejecta due 
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to excessive extinction of redshifted wings of emission lines 
than the blueshifted ones l|Lucv et al.| [l991; Da nziger et al.l 
Il99ll >. Between +300-400d, the observed blueshift in oxy- 
gen may be due to contamination from the Fell mutli- 
plet at 6250A. The reason of the blue-shift of the oxy- 
gen line at early epochs ( < 200 day) is still not clear 
and several hypoth e sis ha ve been done. In a recent work 
iTaubenberger et all (|2009l ) described this as a result of resid- 
ual opacity that remains in the inner ejecta. This seems to be 
the most likely explanation for observed blue-shift of oxygen 
line. Dust formation at an early epoch ~ +300d is also re- 
ported for SN 2004et. For SN 2008gz we estimate a blueshift 
of ~ 250 kms -1 , in [O i] components at epoch of +275d (see 
Fig. [5] and [8]), however due to absence of any other evidence, 
this is not enough to claim dust formation in the SN 2008gz 
ejecta. 



5 PHOTOSPHERIC AND H-ENVELOPE 
VELOCITIES OF EJECTA 

We used multi-parametric resona nce sca ttering code 

SYNOW dBranch. Baron fc Jeffervl |200ll : iBranch et ail 
120021 ; iBaron Il2005h for modelling the spectra of SN 2008gz 
to iterpret spectral features and estimate velocities of lay- 
ers at different epochs. The algorithm works on the as- 
sumptions of spherical symmetry; homologous expansion of 
layers (v ~ r); sharp photosphere producing a black-body 
spectrum and associated at early stages with a shock wave. 
In photospheric phase, the spectral lines are formed by the 
shell above the thick photosphere, but in nebular phase all 
visib le regions are optically thin (|Branch. Baron fe Jeffervl 
l200lf ). Each of these two phases of SN evolution can be ex- 
plained with individual approximations and the modelling 
of observed spectra needed in different synthetic codes. 

Our main aim in modeling the spectral features is to 
estimate the velocities of the layers and that of pseudo- 
photosphere. It was also noted in Branch, Baron and Jeery 
(2001) there is no sharp division between the photosphere 
and the nebular phases. We note the presence of absorption 
components in Iron and Hydrogen lines at the latest spectra, 
which can be explained as a result of the decreasing reso- 
nance scattering mechanism. Although the resonance scat- 
tering codes like SYNOW are not used for describing late 
time spectra (see for ex. Elmhamdi et al., 2006), we use it to 
describe only the absorption parts of line profiles. It is not 
our intention to fit the emission part of the line profiles be- 
cause this procedure needs to make use of other assumptions 
and other codes. 

A preli minary result of SYNO W fit on SN 2008gz is 
reported by iMoskvitin et al.l (|2010l 'l. In Fig. [5] we present 
our model fit for +87d spectrum in detached case, i.e., as- 
suming line forming shells of ionised gases moving ahead of 
the photosphere. Most of the spectral features (particularly 
absorption minima and the continuum) are produced well. 
All the identified spectral features ar e same as ma rked in 
Fig. [5] We tried out undetached cases l |Sonbad l2008) as well 
and also attempted changing density laws (exponential and 
power) and we found that it had very little effect while fitting 
the absorption minima. In order to obtain precise velocity 
measurements of hydrogen layers, we modeled the profiles 
of Ha, H/? and Fen independently (see Fig. 1 1 Of) following 




4000 5000 6000 7000 

Rest Wavelength (A) 

Figure 9. SYNOW modeling of the +87 d spectrum of SN 
2008gz. The line identification is done after considering the SN as 
a spherically expanding fireball and lines are formed in a region 
moving ahead of the photosphere. Optical depth of individual line 
is calculated through "Sobolev approximation" . 
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Figure 10. SYNOW models for absorption features of Hct (left 
panel) and H/3 & Fell (right panel). Blending effect due to Sell, 
Tiii and Ball is also incorporated in the models. Spectral evo- 
lution corresponds (top to bottom) to +87, +115, +140, +171 
+218, +222, +231 and +275. H/3 region is not modeled for +231 
d. 
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Table 6. Velocities of photosphere, Ha and H/3 for different epochs of SN 2008gz evolution. All the parameters are derived from 
SYNOW modelling. 



UT Date 


Phase 


«ph = u(Feii) 


v e (Fe 


n) 


«(Ha) 


i> e (Ha) 






(yy/mm/dd/) 


(days) 


km s — 1 


km s" 


-i 


kms" 


km s 


km s 


-i 


kms 


2008/11/11 


+87 


4200 ± 400 


1000 + 


400 


5500 ± 200 


2100 ± 300 


4300 ± 


300 


1300 ± 100 


2008/12/08 


+115 


3200 ± 400 


1500 + 


500 


4700 ± 300 


1500 ± 200 


3800 + 


600 


900+ 300 


2009/01/03 


+140 


4000 ± 400 


2500 ± 


1100 


5300 ± 400 


1300 ± 300 


4000 + 


400 


1000 ± 200 


2009/02/02 


+171 


3500 ± 300 


1600 + 


800 


5100 ± 200 


1300 ± 200 


3500 + 


300 


1000 ± 200 


2009/03/22 


+218 


3100 ± 500 


1700 + 


900 


4800 ± 400 


1500 ± 500 


3700 ± 


300 


1600 +400 
iouu_ g00 


2009/03/25 


+222 


3100 ± 700 


1200 + 


800 


4900 +700 


1100+ 300 


3100 ± 


700 


2000 ± 1000 


2009/04/03 


+231 


< 5400 






4800+j™ 


1600+ 200 
iouu_ 600 








2009/05/17 


+275 


2000±ioo 


1500 + 


500 


3800 ± 200 


600 ± 200 


2000 ± 


100 


600+ 600 
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Figure 11. Velocity profiles of different elements in the ejecta of 
SN 2008gz 



r ~ exp(—v(r)/v e ), where r is optical depth and v e is the 
c-fold velocity. We also incorporated Till, Sell, and Ball 
ions to model multi-minima absorption features around H/3. 
For Fe II lines, we noticed that velocities of various absorp- 
tion features are similar or have insignificant differences 
and hence we used average d Feu values to represent the 
photospheric velocity (i> p h) (jBranch. Baron fe Jeffervll200ll ; 
lElmhamdi et al.ll2006l ). Estimations of photospheric and en- 
velopes velocities of different layers are given in Table [6] 
Uncertainties in the estimates take into account the noise in 
the spectra. 

We also estimated photospheric and H-envelope veloci- 
ties using IRAF by directly locating the absorption minima 
and the same is shown in Fig. [TT] The velocity for Fe II 4924, 
5018 and 5169 lines, range from ~ 4000 kms -1 at +87d 
to around 3000 kms" 1 at +275d. These values are similar 
to the values estimated above from SYNOW modelling. 
For H/3 and Ha layers our values are consistently higher 
by about 1500 kms -1 at all epochs than that derived from 



SYNOW. This discrepancy may arise due to contamination 
of true absorption minima by the emission component of P- 
Cygni profile of Ha and H/3 and hence it is likely that the 
true blueshift would be overestimated while using absorp- 
tion minima. 

For SN 2008gz the Ha velocity at +87d is ~ 
6800 kms" 1 , while f or SNe 2004et, it is ~ 6000 
kms~ 1 (|Sahu et all I2006J ). H/3 also shows higher expansion 
velocity at comparable epochs. Similarly the photospheric 
velocity at day +50 is ~ 3700 kms" 1 , while it is 4000 
kms^at +87d for SN 2008gz. So, even by considering an 
overestimate of the plateau period about 25 days (i.e. for 
U ~ 95d), the photspheric and H-envelope velocity for SN 
2008gz seem to have comparable or higher values than SN 
2004et. The photosp heric velocity is a good indicator o f the 
explosion energy (see lDessart. Livne fe Waldmadlinol ) and 
hence SN 2008gz has ex plosion energy similar to that of SN 
2004et, ~ 2.3 x 10 51 erg (|Utrobin fe Chugaill2009T ) or higher. 
lUtrobinl (|2007l ) obtains an explosion energy of ~ 1.3 x 10 
erg for SN 1999em, which has comparatively lower expan- 
sion velocity than to SN 2004et at similar epochs. 



6 DISTANCE AND EXTINCTION OF SN 
2008GZ 

Spectrum of the nuclear region of the host galaxy taken 
on +170d (see Fig. Qjffl) was used to estimate czheiio, t ne 
heliocentric velocity. Employing five nebular emission lines 
and prominent absorption features and using the 5577 
A skyline as a reference wavelength, we obtain c^eUo of 
1891 ± 14 kms -1 This is in agreement with 18 other ra- 
dio and optical measurements of czhciio of NGC 3672 (in 
range of 1400 to 2000 kms" 1 ) listed in HyperLEDA0- The 
Combined measurement gives a mean czhciio of ~ 1864 ± 19 
kms" 1 and it corresponds to a corrected (Local Group infall 
into Virgo) distance of ~ 25.65 ± 2.93 Mpc0- NELO lists 
four distance measurements based on H I Tully-Fisher rela- 
tion, with a mean of 25.25 + 4.0 Mpc, which is in agreement 



8 Fig. I12l is available only in electronic form. 

9 http://leda.univ-lyonl.fr/ 

10 The cosmological model with Ho = 70 km.s _1 .Mpc _1 ,f2 m = 
0.3 and f2jv = 0.7 is assumed throughout the paper and the un- 
certaint y corresponds to a local cosmi c thermal velocity of 208 
kms- HTerrv, Paturel fe Ekholmll2002h . 

11 http:/ /nedwww. ipac.caltech.edu/ 
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Figure 13. Rest wavelength spectra of SN 2008gz (+87 d/lOA 
and +115 d /6A) and the centre of host galaxy (+170 d /6A). 
Na I D absorption due to interstellar matter of host galaxy (~ 
5892A) and the Milky Way (~ 5854A) is indicated. 



Table 7. Equivalent width measurement of Na I D absorptions in 
the spectra of SN 2008gz and the host galaxy. Last row provides 
the uncertainty weighted EW of Na I D absorption in the direction 
of SN 2008gz due to Milky Way and the host. 



UT Date 


Phase 


EW (MW) 


EW(host) 


(yy/mm/dd/) 


(days) 


A 


A 


2008/11/11 


+87 


1.39 + 0.34 


0.23 ± 0.28 


2008/12/08 


+115 


1.21 + 0.69 


0.28 ± 0.50 


2009/02/02 


+ 171 


1.32 ± 1.29 




Weighted EW 




1.29 + 0.29 


0.23 ± 0.24 



with the above kinematic estimate, we therefore, adopt un- 
certainty weighted distance of 25.5 + 2.4 Mpc for SN 2008gz. 

The Galactic reddening in the direction of SN 2008gz 
as derived from the 100 /im all sky dust extinction map of 
ISchleeel et "all ljl998h is estimated as E(B - V) = 0.041 ± 
0.004 mag. Additionally, we could also determine redden- 
ing in the direction of SN 2008gz from equivalent widths of 
Na I D absorption lines present in the spectra of SN 2008gz 
(+87d and +115d) and the centre of the host galaxy (see 
Fig. \\^. Di (5889.95A) and D 2 (5895.92A) component of 
Na I D is not resolved in +87d spectra and it is seen as nar- 
row absorption features overlaid on the broad P-Cygni emis- 
sion wings of NalD due to SN. In rest wavelength plot the 
host galaxy contribution is seen at ~ — 80kms _1 , while the 
Galactic contribution is at — 1630 km s -1 . In +115d spec- 
trum, both components of NalD are resolved. Intriguingly, 
the Galactic component appear to split into two - a stronger 
component at ~ - 1896 km s" 1 (due to Milky Way ISM) 
and a weaker component at ~ — 1021kms _1 , possibly due 



to inter-galactic medium. In the +170d spectra, the Galac- 
tic NalD absorption appear as a single component. Esti- 
mated total NalD equivalent widths (EW) are reported in 
Table [7] Quoted errors in EW are photon-noise dominated 
RMS uncertainty derived following I Vollmann fc Eversbergl 
(2006, see their Eq. 6). It is seen that the EW contribution 
due to host (0.23 + 0.24A) is smaller in comparison with the 
total Galactic contribution (1.29 ± 0.29A). 

It is known that the EWs of interstellar absorption 
bands is well correlated with the reddening E (B — V) esti- 
mate d from the tail o f SNIa colour curves (|Barbon et alj 
Il99d : iRichmond et al.l Il994|) and by employing empiri- 
cal re lations established by iTuratto. Benetti fc Cappellarol 
(|2003h . E(B — V) = -0.01 + 0.16EW (where EW is in A0, 
we obtain Galactic E(B — V) contribution as 0.20 ± 0.05 
mag and host galaxy as 0.03 ± 0.04 mag. The Galactic 
E(B — V) derived in this way is larger than that derived 
from Schlegel map. Considering the normal extinction law 
(Rv = 3.1) and Schlegel value of E(B - V) = 0.041 for the 
Milky Way; the EW (host/Galactic) ratio would suggest a 
slightly lower reddening in the host {E(B — V) < 0.01 mag). 
However the host galaxy value of 0.03 mag may not be ruled 
out in case of a different dust to gas ratio for the host. 

For this work, we will adopt a conservative value of 
E(B - V)=0.07 ± 0.04 mag, obtained by adding Galac- 
tic (Schlegel) and host galaxy (NalD) contribution. This 
corresponds to visual extinction (Ay) of 0.21 ± 0.12 by 
assuming ratio of t otal-to-selective extinction Rv = 3.1 
ICardelli et al.ll 19891 ). 



7 TEMPORAL EVOLUTION OF COLOUR 
AND BOLOMETRIC LUMINOSITY 

Fig. [14] shows reddening corrected colour evolution of SN 
2008gz. For comparison, we also show reddening corr ected 
colour c urve of SN 1987A dsiintzeff fc Bouchet|ll990l), SN 
1999e m l|Elmhamdi et al.ll2003aT ) and SN 2004et (|sahu et all 
I2006I ). for E(B-V) of 0.15, 0.10, and 0.41 mag and explosion 
epochs of JD 2446849.82, JD 2451480.5 and JD 2453275.5 
respectively. Though the colour curves of SN 2008gz have 
large scatter but overall nature of its temporal evolution is 
prominent. SN 2008gz follows general trend of colour evo- 
lution i.e. a steep and rapid decrease from blue (high tem- 
perature) to red (low temperatures) colours similar to SNe 
2004et and 1999em. (B — V)o becomes redder from 1.2 mag 
at +87d to about 1.7 mag at +140d and it follows trend 
similar to SNe 2004et and 1999em. The overall trend in 
colour evolution of SN 2008gz between end of plateau and 
middle of the tail is similar to that of peculiar type II SN 
1987A. In nebular phase, (B — V)o turns blue rather rapidly, 
— 1 mag (lOOd)" 1 and this arises due to suspected flattening 
in B light and a shallow decay in V, R and I bands. The 
(V — R)o colour is found to be consistently bluer than for 
other IIP SNe. For (V — R)o and (V — I)o, the increment is 
quite shallow during the transition from plateau to nebular 

12 In ITuratto. Benetti fc Cappellarol d2003ri there are two rela- 
tions — one is with low slope and other with a high slope. In 
this work we have considered the lower slope, because it is well 
sampled and match es with the previous works in this direction 
llBarbon et alJll99Ch . 
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Figure 14. Temporal variation of colour of SN 2008gz. Also 
shown are the other core-collapse supernovae, SN 1987A, SN 
1999em and SN 2004et. 
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Figure 16. Bolomctric light curve of SN 2008gz. For comparison, 
we also show light curves of SN 1987A, SN 1997D, SN 1999em 
and SN 2004et. 
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Figure 15. Flux contribution in percentage, in BVRI bands of 
SN 2008gz along with a comparison to SN 1987A 

phase. This is in contrary to the colour evolution of low lumi- 
nosity type II-peculiar SNe 1997D and 1999 eu, where a steep 
rise and excess in colou r has been noticed (|Pastorelk> et al.l 
l2004lMisra et alj|2007ft . 

Bolometric luminosity is essential to estimate total opti- 
cal radiant energy in the explosion and also at the tail phase 
it is a good estimator of radioactive 56 Ni synthesized in the 



explosion. To a good approximation, the integration of ex- 
tinction corrected flux in UBVRI at a given epoch gives a 
meaningful estimate of bolometric luminosity. The extinc- 
tion corrected BVRI ma gnitudes were converted into f luxes 
using zeropoints given bv lBessell. Castelli fc Plezl (|l998h and 
the total BVRI flux is obtained by interpolating and inte- 
grating fluxes between 0.4 to 0.85 /J,m. Fig. ll5l shows percent- 
age flux contribution in different passbands and the overall 
trend is found to be similar to that of SN 1987A. From 
plateau (+87d) to nebular (+140d) phase, the flux contri- 
bution declines from ~ 15 to 10% at B, from ~ 25 to 20% at 
V, while it ascends from ~ 36 to 42% at R, from ~ 22 to 30% 
at I. In nebular phase (until +270d), the flux contributions 
remain constant. For nearby SNe 1987A and 2004et, large 
(~ 40%) flux contribution in U and B bands are observed 
at initial epochs, which reduces to a few percent by +60d. 
Though, towrads later epochs, when SN ejecta becomes op- 
tically thin, a little enha ncement (about 5 %) in U and B 
bands is also noticeable l|Misra et alj 120071. Therefore, for 
SN 2008gz, we have constructed UBVRI bolometric light 
curve after making a constant (5%) contribution from U 
band over the period of our observation. No correction for 
flux contributions in the ultraviolet and near-infrared region 
were made as they become significant respectively in early 
and late phases of the light curve evolution. 

Fig. [16] shows the UBVRI bolometric nature of SN 
2008gz along with other type II events. Different behavior 
of these events is clearly evident and it provides constraint 
on the synthesized radioactive 56 Ni as well as explosion en- 
ergy of SNe. Tail luminosity of SN 2008gz is similar to SN 
2004et, while the plateau luminosity is about 0.2 dex fainter. 
Explosion parameters for SN 2008gz are estimated in next 
sections. 
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8 PHYSICAL PARAMETERS 

8.1 Amount of ejected radioactive Nickel 

The nebular phase light curve of type II SNe is mainly gov- 
erned by the radioactive decay of 56 Ni to 56 Co to 56 Fe hav- 
ing half- life of 6.1 and 111.26 days respectively and hence 
the tail luminosity is directly proportional to the amount of 
56 Ni synthesize d by explosive burning of S i and O during 
shock breakout (|Arnettlll980l ; lArnett 1119961 ). 

By using tail luminosity, 56 Ni mass can be derived fol- 
lowing the method described by iHamuvl |2003l) applied un- 
der the assumption that all the 7-rays emitted during the ra- 
dioactive decay make the ejecta thermalised. For SN 2008gz 
using the V band magnitude at +200d, corrected for extinc- 
tion (Ay = 0.21 ± 0.12 mag; Et , a bolometric correction of 
0.26 ±0.06 mag (|Hamuvl l2001') during nebular phase, and a 
distance modulus of 32.03±0.21, we derive tail luminosity of 
1.51 ± 0.29 x 10 41 ergs -1 and this, for the plateau duration 
of 115 days, results in Ni mass M Ni = 0.067 ± 0.012M Q . 

We can estimate the mass of 56 Ni by a direct com- 
parison with that of SN 1987A for which it is accurately 
determined as O.O75M0. We note that the temporal evolu- 
tion of flux contribution in different bands for SN 2008gz is 
comparable with that for SN 1987A and consequently 
we can safely assume that at a comparable epoch the ra- 
tio of their luminosities is equal to the ratio of synthesized 
56 Ni mass. SN 2008gz attains the deep nebular phase be- 
yond 150 days after the burst. The last observation was also 
performed nearly at +275 day. Hence, the mean ratio be- 
tween tail UBVRI bolometric luminosity of SN 2008gz (~ 
150-275) and that of SN 1987A is about 0.539. This im- 
plies that for SN 2008gz ejected 56 Ni mass [0.539x0.075] w 
0.041M Q . 

By taking a sample of ten IIP SNe, lElmhamdi et al.l 
( 2003b) show that the steepness of V-band light curve slope 
(defined as S=dmy/dt) at the inflection time (t;) is anti- 
correlated with 56 Ni mass (logA'/ N i = -6.2295S - 0.8147). 
For SN 2008gz we have well sampled transition phase and 
get a value of S = 0.075 ± 0.036, mag d" 1 (see Fig. ITfFh 
which corresponds to Mjvi = 0.052 ± O.OIMq. Considering 
the uncertain values of extinction to wards SN 2008gz , we 
note that for SNe 2004A and 2003gd, iHendrv et~aD (|2006l ) 
finds that the lElmhamdi et al.l l|2003rj ) scheme gives some- 
what lower values. Further, based on plateau luminosity, a 
lin ear correlation, logM N i = -0A38M V (U -35) -8.46 found 
by lElmhamdi et al.l l|2003bh provides a value of 56 Ni mass 
as O.O51M for M v of -16.37 ±0.24 at ±87d (U - 28). Tak- 
ing average of above four estimates, we get the amount of 
produced 56 Ni mass = 0.05 ± O.O1M . For estimation of 
other physical quantities we assume that above amount of 
56 Ni was produced by SN 2008gz. Though we would like 
to state that adopting smaller plateau duration ( t]2.2[) will 
reduce the amount of ejected radioactive 56 Ni which will 
further propagate in determination of progenitor properties 

(S3}. 



8.2 Environment of the progenitor 

Constraints on the nature of progenitors of core-collapse 
SNe are derived from the studies of environments in which 
they occur. For example, by correlating the position of 
explosion sites with that of the sites of recent star for- 
ma tion as traced by Ha l i ne emission, it is found (see 
e.g.lJames fe Andersonll2006l ; lKellv. Krishner fc Pahrell2008l ; 
lAnderson fe Jamesll2008l ) that core- collapse events are excel- 
lent tracer of star formation, and type Ib/c are more likely to 
be associated with regions of hig h surface brightness o r high 
Ha emission than type II SNe. lAnderson fe Jame s (2009) 
found that type IIP events are likely to be more centrally 
concentrated than other II-sub types. SN 2008gz occurred in 
spiral arms of the host galaxy at a deprojected galactocentric 
distance of 2.8 kpc (within half light radius) and the oxy- 
gen abundance ([0/H] = 12±logjVb/./VH) of the galactic ISM 
at the position of SN is estima ted as 8.6 (derived from the 
O/H-M b relation proposed by iPilvugin. Vflchez fc Continil 
2004) ), which is clos e to the solar abundance [O/H] of 8.65 



Asplund et all 12009). 



13 Fie. 1171 is only available in electronic form. 



In order to further investigate the level of Ha emission 
level, we observed NGC 3672 in narrow-band Ha-line (A c 
= 655lA) and Ha-red (A c = 6650A) having FWHM of 83A 
and 79A respectively. We used 2k CCD camera mounted 
with lm ST, Nainital on 18 February 2010 (±560 d). A to- 
tal exposure of 1 hr in each Ha-line and Ha-red filters were 
taken along with several bias and sky flats. Raw images were 
corrected for bias and flats using IRAF. FWHM of stellar 
PSF (seeing) varied from l'.'8 to 2'.'2. Images were combined 
to improve signal-to-noise ratio and the continuum subtrac- 
tion was done using ISIS. Our narrow-band filter-set were 
not customized for extra-galactic work, and we expected 
Ha-red to contain emission line fluxes, as at the redshifted 
wavelength (~ 4lA at Ha) of NGC 3672, the Ha-red fil- 
ter had a transmission of 25%, 40% and 80% respectively 
for Nll6548A , Ha and Nll6589A. We could verify this by 
subtracting Ha-line frame from the broad-band R and V 
frames of SN taken on 14 February 2010, which gave no 
residual, while Ha-red frame showed residuals. Fig. ll8l shows 
the continuum subtracted image of SN 2008gz showing con- 
tributions from Ha±Nll emissions. A close up view of SN 
location is also shown. The sky is at level of while the 
peak flux is around 22 counts. Five prominent regions of Ha 
emission having peak counts above 15 are clearly apparent. 
SN 2008gz position is at a level of 8 counts, and hence it 
belongs to a low-luminosity H 11 regions. This is also evident 
from the early epoch spectra, in which an emission peak (not 
so prominent) of zero velocity is seen. 



8.3 Properties of progenitor star 

Accurate estimates of explosion parameters require detailed 
hydrodynamical modeling of the optical light curves, though 
the analytical relations (based on a few well modeled IIP 
events) correlating the physical parameters explosion energy, 
pre-SN radius and total ejected mass on the one hand and 
the observable quantities, plateau duration, mid-plateau V- 
band magnitude ((My) mp ) and mid-plateau photospheric 
velo city (fm P ) on the other hand are proposed to exit (see 
e.g. iPopovl 1 19931 ; iLitvinova fc Nadvozhinl Il985l ; iNadvozhinl 
2003). For SN 2008gz , we can only provide an approximation 
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Figure 18. Continuum subtracted Ha image of NGC 3672. In left panel (~ 3' X 3') five bright Hn regions are visible. The SN position 
is marked with a circle. In right panel (~ 1' X 1') zoomed in image of SN is shown. SN location is marked with a dot. 



of these observables. The estimates of (Mv)m P and v mp can 
be made by using the 56 Ni mass estimate and by employing 
following empirical relations d erived (by m inimizing least 
squares) from the data given in lHamu v (2003, see their Fig. 
3 and 4). 

log(M Ni ) = -0.385 x (Afv)mp - 7.749 
log(M Ni ) = 2.771 x log(« mp ) - 11.425 



For M Ni = 0.05 ± 0.01 M , we derive (M v ) r _ 
— 16.7^3 mag and « mp = 4503 



+306 
348 



km 



This 



value of absolute magnitude is consistent with those ob- 
tained from the first photometric V point (£JHJ). These es- 
timates, along with a plateau duration of 115 days provides 
jLitvinova fc Nadvozhin]|l985t ) burst energy ~ 2.5t° S 7 x 10 51 
erg, ejected mass ~ 34tg° Mq and pre-SN radius ~ 167^gi 6 
R.0. The explosion energy derived in this way is consistent 
with the one expected from photospheric velocity (see fJSJ) , 
however, the ejecta mass is larger in comparison to the typ- 
ical progenitor mass range (8 .5-16.5Mn) estimat e derived 
from pre-explosion SN images (|Smartt et alj|2009al ) and this 
may arise due to uncertainty in above m easured parameters. 

iDessart. Livne fc Waldmanl (|2010f ) demonstrate that by 
employing the explosion energy estimate, the observed line 
width of [O i] and artificially generated radiation hydrody- 
namic simulations of core-collapse SNe, it is possible to put 
an upper limit on the main-sequence mass of the progenitor. 
For SN 2008gz, the Oxygen ejecta velocity of ~ 1250 km s _1 
(HWHM of [Ol] profile, see © (which is slightly higher 
than ~ 1000 km s" 1 observed for SNe 2004et and 1999em at 
around +330 d) and assuming an explosion energy of 3 x 10 51 
erg s _1 provide a main-sequence mass of 15 (12) respectively 
for non-rotating (rotating) pre-SN models while assuming an 
ejecta velocity of 1500kms _1 for [Ol] gives an upper limit 
of 17(13) Mq, which is consistent with t hat derived using 
pre-explosion images jSmartt et alj|2009ah . 



-10 





I ' 


I ' 


I 

. 2008gz 








□ 2004et 








o 2004 dj 








4 2004A 








□ 2003gd 


* 






» 1999em " 
o 1999gi 








o^j 1998A ~ 
o 1997D 




o <%, 1 

°® ° 




, 1 992] 1 

« 1990E 

19B7A - 

— . "Co -> "Fe 




°© 

o 


% 


D 






o 

o 


□ 




1 


o 

I 


1 





100 200 

Days 


300 400 
since explosion 


500 



Figure 19. Comparison of absolute V-band light curve of SN 
2008gz with other type IIP SNe, like SN 2004et, SN 2004dj, SN 
2004A, SN 2003gd, SN 1999em, SN 1999gi, SN 1998A, SN 1997D, 
SN 1992H, SN 1990E and SN 1987A. The magnitudes have been 
corrected for distance and readdening. 



9 COMPARISON WITH OTHER CORE 
COLLAPSE SUPERNOVAE 

A detailed investigation of SN 2008gz indicates that it is 
a normal type IIP event showing photometric and spec- 
troscopic evolutions similar to archetypal SNe 2004et and 
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1999em. SN 2008gz occurred in a highly inclined (Oi nc = 
56.2°) host galaxy, within deprojected galactocentric radius 
of 0.27r25 (Table[T]) implying solar metallicity region similar 
to SNe 1999em and 2005cs, though its explosion and other 
properties were found to be similar to SN 2004et, which oc- 
curred in the outskirt of its host galaxy. Our narrow-band 
Ha photometry indicates that SN 2008gz was associated 
with a star forming low luminosity H n region of the galaxy. 
Thus, the metallicty appear to have little effect on the ex- 
plosion properties of core-collapse SNe. 

In Fig. [19] we show absolute V-band light curve of well 
studied core-colla pse events collected from literature (see 
iMisra et al.l 120071. for references) . The plateau luminosity 
(M v ~ -16.60 of SN 2008gz is at the similar level as 
of peculiar SNe 1987A, 1998A and it lies in between to the 
brighter end of IIP SNe 1990E, 1992H, 2004et (~ -17 mag) 
and the fainter end SN 1997D (~ —15 mag). This shows that 
SN 2008gz is a normal event, both in energetics and nickel 
production. These values of plateau luminosity for IIP SNe 
are lower than —17.6 ± 0.6 mag w hich was predi cted us- 
ing theoretical models calculated by Irloflichl (|200lh for IIP 
SNe based on wide range of parameters (explosion energy, 
metallicity, mass loss of progenitor). 

SN 2008gz showed rarely observed 1.5 mag drop at V 
from plateau to nebular phase and it had tail luminosity 
comparable (or higher) to SNe 2004et resulting in the syn- 
thesised 56 Ni mass in range 0.05 — O.IM©. Large tail lumi- 
nosity of SN 1998A indicates that the thermalization pro- 
cess for 1998A was more efficient than SN 2008gz. Colour 
evolution of SN 2008gz has a similar trend like normal IIP 
and peculiar SNe 1987A and 1998A. Expansion velocity 
of the SN 2008gz ejecta was comparable to SN 2004et or 
higher, implying explosion energy of ~ 2 x 10 erg. Our 
calculation for progenitor mass ( H8.3[) of SN 2008gz favours 
the mass range of 13 - 18 M (i.e., SNe 1999em, 1999gi, 
2004dj, 2004et) than to the lower mass range of 8 — 18 
Mq (i.e., SNe 1997D, 2004A, 2005cs). This progenitor mass 
grou ping is also favoured on the basis of radio luminos- 
ity (|Chevalier et al.ll2006l ). As an alternative scenario, SN 
2008gz can be characterized as a peculiar event. Lack of data 
in first 1-2 month s after the explosion and spectral similarity 
with SN 1998A |Benetti et al J l2008h also indicates toward 
the possibility that SN 2008gz is a peculiar type II event. 
Figure [20] shows a comparison of SN 2008gz spectrum with 
SN 2004et and SN 1998A at comparable epoch. There are 
many similarities between spectral features of SN 2008gz 
and SN 1998A. On the other hand, low rates of peculiar 
type II SN, similarity in late light curve of SN 2008gz with 
type IIP events and similarity of some particular spectral 
features, like high velocity Ha line (Q of this event with 
that of SNe 2004et and 1999em also indicates that SN 2008gz 
would be a normal type IIP event. 



14 We estimate the value of mid plateau My as — 16.6±0.2 mag, 
by considering an average decline rate of 0.006 mag d —1 during 
the plateau phase of SNe 1999em and 2004et, we found that the 
mid-plateau magnitude of SN 2008gz was ~ 0.2 mag brighter 
than that determined at +87d (SUTTJ. We assume Ay of 0.214 
mag. 
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Figure 20. Comparison of spectral features of SN 2008gz along 
with SN 2004et and SN 1998A 



10 SUMMARY 

We present BVRI photometric and low resolution spectro- 
scopic observations of a supernova event SN 2008gz which 
occurred in an spiral arm and within half-light radius of 
a nearby (~25 Mpc) galaxy NGC 3672. As the event was 
buried the galaxy light, we used template subtract technique 
to estimate the apparent magnitude of the event. We mon- 
itored the event for a period of ~ 200 days. We summarize 
our results as follows. 

(i) Photometric and spectral nature of the event is found 
similar to normal type IIP SNe 2004et and 1999em. The 
event was discovered by about 82 days after the burst and 
it has a plateau phase lasted 115 ± 5 days. We monitored 
the SN evolution from +87d to +275d. 

(ii) We estimate photospheric and H-envelop velocity by 
using both direct measurements of the absorption minimna 
of H I and Fe n lines and SYNOW modelling of the spectra. 
Here both values agree well within uncertainties. We esti- 
mate the photospheric velocity of ~4000 kms _1 at +87d, 
which is higher than that observed for well studied SN 
2004et at similar epochs indicating explosion energy com- 
parable to or higher than 2004et. Similar trend was also 
seen for the expansion velocity of H-envelopes. 

(hi) Using pre-SN models of iDessart. Livne fc Waldmanl 
(2010) and also by comparing explosion energies (derived by 
using hydrodynamical models) of well studied IIP SNe, we 
find that the SN 2008gz had explosion energy of 2 - 3 x 10 51 
erg.s -1 . This estimation, coupled with the observed width 
of the forbidden [O i] line gives an upper limit for the main- 
sequence progenitor mass of 17Mq. 

(iv) SN2008gz exhibits rarely observed drop of 1.5 mag 
within 30 days in V-band from plateau to nebular phase, 
this is higher than the typically observed fall of 2-3 mag 
in normal IIP SNe. Adopting Av = 0.21 mag, we could 
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estimate mass of 56 Ni synthesized during the explosion as 
0.05 ±0.01M Q . 

(v) Our Ha observation taken about 560 days after the 
explosion indicates that the event took place in a low lumi- 
nosity star forming arms, very close to a Hn region. The 
emission kink of this H II region is also seen in Ha line near 
zero velocity of the Doppler corrected spectra of SN. 
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Table 2. Journal of photometric observation of SN 2008gz. 
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a With reference to the time of explosion JD 2454694.0. 

6 LT : 0.6 m Lowell Telescope, Perth Observatory, UK; ST : 1 

m Sampurnanand Telescope, ARIES, India; IGO : 2 m IUCAA 

Girawali Observatory, India; NTT : 3.6 m New Technology 

Telescope, ESO, Chile 

c FWHM of stellar PSF at V band 

d Ellipticity is estimated using IMEXAMINE task of IRAF 
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Figure 3. Comparison of ISIS derived magnitudes and ours for BVRI. 
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Figure 8. Two component Gaussian fit of [O i] 6300, 6364 A emission lines. Measurements of their relative intensities quantifies, whether 
the corresponding line emitting region is optically thick or thin. 



22 Roy et al. 



1.5 



x 

. — i 

T) 

0) 
N 



o 
Z 



0.5 




4000 



5000 6000 
Observed Wavelength (A) 



Figure 12. Spectrum of nucleus of the host galaxy NGC 3672 taken with 2-m IGO. The spectrum shows Ha, Nil (6548 A and 6583A), 
and Sn(6717Aand 6730A) lines in emission (similar to Sc spiral galaxies), which can arise from a gaseous component heated by AGN, 
post-AGB stars, shocks or cooling flows. The Call(K and H), H7, H/3, MglT, Fel(5270A), NalD due to Milky Way and the host are 
seen in absorption. 




Figure 17. Steepness parameter estimation for SN 2008gz. 



